Finite temperature effects in trapped Fermi gases with population imbalance 
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We study the finite temperature T behavior of trapped Fermi gases as they undergo BCS-Bose Einstein con- 
densation (BEC) crossover, in the presence of a population imbalance. Our results, in qualitative agreement with 
recent experiments, show how the superfiuid phase transition is directly reflected in the particle density profiles. 
We demonstrate that at T ^ and in the near-BEC and unitary regimes, the polarization is excluded from the 
superfiuid core. Nevertheless a substantial polarization fraction is carried by a normal region of the trap having 
strong pair correlations, which we associate with noncondensed pairs or the "pseudogap phase". 



PACS numbers: 03.75.Hh, 03.75.Ss, 74.20.-z 

Recent work Ql 0, 01 on trapped atomic Fermi gases with 
population imbalance has become particularly exciting. With 
the application of a magnetic field, these systems exhibit a 
continuous evolution \4, 5, 6]fromBCS to Bose-Einstein con- 
densation (BEC). Not only are these gases possible prototypes 
for condensed matter systems |7, 8] in the presence of a mag- 
netic field - Zeeman coupling, but they may also be proto- 
types for particle and nuclear physics systems l9l. [l(jll . These 
pioneering experiments have been done so far by two experi- 
mental groups 1 1, 2]. And what differences are present appear 
to lie more in the interpretation than in the actual data. 

There are a number of key experimental observations which 
we now list, (i) Both groups have observed that the trap pro- 
files are characterized by a central core of (at most) weakly 
polarized superfiuid, surrounded by a normal region where 
the bulk of the polarization is contained, (ii) The normal 
region appears to consist of overlapping clouds of both spin 
states ("normal mixture"), followed at the edge of the cloud 
by a region consisting only of the majority component. There 
is not complete agreement |2J, |3[ on whether the normal- 
superfiuid boundary is sharp which would correspond to some 
form of phase separation. 

These population imbalance experiments have been done 
01 IH in conjunction with other measurements (vortex exci- 
tations and magnetic field sweeps) which establish the pres- 
ence and the location for superfiuid condensation, (iii) Even 
more recently |3] it has been demonstrated that the superfiuid 
phase transition at T c can be directly reflected in changes in 
the shape of the clouds, (iv) Important for the present pur- 
poses is the fact that |3] there are strong interaction effects 
within the normal region of the cloud. 

The goal of the present paper is to address the four points 
[(i)-(iv)] listed above through a finite temperature theory of 
BCS-BEC crossover in the presence of population imbalance 
within a trap. A related study of the homogeneous system 
was presented earlier ITTIl . What is unique to our work is the 
capability of separating in a natural way the condensed from 
noncondensed pair contributions to the trap profile. The diffi- 
culty in making this separation lies in the fact that (except at 
T = 0) the presence of a fermionic excitation gap, is not a sig- 
nature of phase coherent superconductivity. We also present 
calculations of T c in a trap and show how the general shape of 
the profile changes below and above T c , unlike what is found 
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experimentally and theoretically fl2^ in the absence of polar- 
ization. In a related fashion, we examine the noncondensed 
pair states in the trap and determine to what extent they differ 
from a free gas mixture of the two spin states. 

Our principal findings are at general T ^ and for the uni- 
tary and near-BEC regimes, (a) the superfiuid core seems to be 
robustly maintained at nearly zero polarization, (b) The mixed 
normal region, carries a significant fraction of the polarization 
within a non-superfluid state having strong pair correlations. 
Indeed, experiments suggest 1 3| that "even in the normal state, 
strong interactions significantly deform the density profile of 
the majority spin component". Here we interpret these cor- 
relations as noncondensed pairs which have no counterpart at 
T = and which are associated with an excitation gap ("pseu- 
dogap") in the fermionic spectrum. Finally, (c) in the course 
of making contact with points (i)-(iv) listed above we show 
good qualitative agreement with experiment. 

Because we restrict our attention to condensates (and their 
pseudogap phase counterparts) with zero momentum (qo = 0) 
pairing, we do not explore those regimes of the phase dia- 
gram corresponding to the lowest temperatures, and highest 
polarizations. Recent very nice theoretical work based on the 
Bogoliubov-de Gennes (BdG) approach has shown 

that in the ground state at unitarity the qo ^ 0, Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) state 1 15] must be incorporated. 
Importantly, the polarization in this state at T — appears 
at the edge but within the condensate 11411 . Fortunately, the 
present work provides a good indication of where the FFLO 
phase will enter, since it occurs when the qo = phase is 
found to be unstable 

dm- 

With the finite T required by 
experiments, it is quite possible that the FFLO phase presents 
itself, primarily in the form of noncondensed pairs, which lie 
on the perimeter of the condensate. 

The value of the present work is derived from the fact that 
a central theme in the experimental literature involves distin- 
guishing the condensate from the normal regions of the trap. 
Whether there is phase separation or not Ql 0] the precise 
nature of the normal (N) and superfiuid (S) phases are all of 
great interest. In this way one needs a theory which distin- 
guishes N from S at finite temperatures, where the excitation 
gap is no longer a signature of superfluidity. Previous theoret- 
ical approaches, based on the Bd G fl3l[l4ll and local density 
approximation (LDA) fl7t d, Q^, l20ll schemes have empha- 
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Figure 1: (Color online) Spatial distribution of the excitation gap A(r) (black) and order parameter A ac (r) (red, main figures) and density 
n CT (r) (insets) at 1/kpa = 1.5 for the majority (black) and minority (red) fermions at different polarizations (p = 0.15, 0.35, and 0.7 from 
left to right) above (upper row) and below (lower row) T c . Here T c /Tf ~ 0.36 ,0.35, and 0.31, respectively. The density difference Sn(r) is 
shown in blue in the main figures, sharing the same vertical axis as A(r). The temperatures for the upper row are T/Tf = 0.4, 0.4, and 0.35, 
respectively, and for the lower row T/Tf = 0.1. Shown on the far right is the difference in column density, Sri2d(r), for the three polarizations 
above (red) and below (black) T c . Here Ef = h 2 k F /2m = UbTf is given by the Fermi energy for an unpolarized, noninteracting Fermi gas 
with the same total number N at T = 0, and Rtf = \J1Ef /mu 2 is the Thomas-Fermi radius. The units for n and 8n2d are k% and k F , 
respectively. 



sized T = 0, albeit without reaching any clear consensus. The 
inclusion of finite T for the LDA case has been introduced 
within the same formalism we use here @], but without sepa- 
rating the condensed and noncondensed pair contributions. In 
addition the application of BdG to T ^ is viewed as prob- 
lematic because it does not incorporate noncondensed pairs 
At the same time, it should be stressed that this 
BdG approach 1131 1 1411 is most likely the appropriate way to 
get a full picture of the T = phase. 

The formalism used in this paper was outlined earlier full . 
Here we incorporate trap effects by use of the LDA. We adopt 
a one-channel approach since the 6 Li resonances studied thus 
far are broad and consider a Fermi gas of two spin species with 
kinetic energy = h 2 k 2 /2m subject to an attractive contact 
potential (U < 0) between the different spin states. We define 
5n = — > 0, where n — n-[ + ri| is the total atomic 
density. Importantly, we include 12311 noncondensed pairs at 
general T. The T-matrix or noncondensed pair propagator, 
is t(Q) = U/[l + Ux(Q)], where x(Q) i s me pair suscep- 
tibility discussed earlier which depends self consistently on 
the fermionic excitation gap A. The presence of pairing cor- 
relations means that A 2 contains two additive contributions 
from the condensed (A 2 C ) and noncondensed pairs (A 2 ff ). In 
the superfluid phase, we have 1 + U x(0) = 0, equivalent to 
Impair = 0, the BEC condition of the pairs. As a consequence, 
the equations become simpler below T c and we may expand 
the T-matrix to arrive at a characteristic frequency J7 q which 
characterizes the dispersion of the noncondensed pairs. 

We now summarize the self-consistent equations JOHHl, 
in the presence of a spherical trap, treated at the level of LDA 
with trap potential V(r) — \mio 2 r 2 . T c is defined as the 



highest temperature at which the self-consistent equations are 
satisfied precisely at the center. At a temperature T < T c the 
superfluid region extends to a finite radius R sc . The particles 
outside this radius are in a normal state, with or without a 
pseudogap. 

The generalized local gap equation is given by 
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where n pa i r (r) = in the superfluid region r < R sc , and 
must be solved for self-consistently at larger radii. The quan- 
tity Z is the inverse residue of the T-matrix 1 1 lj. For conve- 
nience we write f(x) = [f(x + h) + f(x — h)]/2 where f(x) 
is the Fermi distribution function. Here we set H = 1, and 
use m/Aira = 1/U + X)k(^ ek ) _1 t0 regularize the contact 
potential, where a is the two-body s-wave scattering length. 
The dispersion is given by E k = — /i(r)] 2 + A 2 , with 

fi(r) = (fi^+fj,i)/2 — V(r). We also define the r-independent 
parameter h = (/ij — /ix)/2. Since Sn > 0, we always have 
h > 0. More generally, [i a is the chemical potential for spin 
a at the trap center. 

The local pseudogap contribution (present only at T ^ 0) 
to A 2 (T) = A 2 SC (T) + A 2 (T) is given by 
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where b(x) is the usual Bose distribution function. The den- 
sity of particles at radius r can be written as 



(3) 
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which depends on the coherence factors u^.,v^ = (1 ± 
6c/^k)/2 with Ck = £k - /■*(?")> and = —h + Ek ,and 
i?ki = h + £k- The total number of particles and the polar- 
ization are respectively given by 

N a (r) = Jd 3 rn a (r), N = + N u (4) 
p = (N^+NJ/N. (5) 

Our calculations proceed by numerically solving the self- 
consistent equations. Here we use the N and p as input; these 
are the control parameters in experiment. 

Figure n shows the behavior of the various gap parameters 
and the majority and minority spin components as a function 
of radius in the trap, for the case of a near-BEC system with 
1/kpa — 1.5. The upper panels are for the normal phase 
and the lower panels are in the superfluid state. We present 
results for three different polarizations and focus first on the 
lower panels where there are two distinct components to the 
gap A sc and A pg . The two gap functions, A sc and the to- 
tal gap A, are plotted vs r along with the difference density 
for up and down spins, or alternatively, the polarization. We 
overlay these plots into order to show clearly what are the con- 
tributions to the polarization from the condensate (I), where 
A sc 7^ 0, the correlated, but normal mixed region (II), where 
A sc = 0, but A 7^ 0, and non-interacting Fermi gas(s) regime 
(III), where A = 0. 

It can be seen that there is very little polarization present 
in the condensate (I) which appears below T c , as has been 
inferred experimentally iQJ y|, y| • Rather the bulk of the po- 
larization is present in the correlated, but normal region (II) in 
which there is a finite excitation gap A, but vanishing A sc . In 
region III at even larger radii, A is essentially zero and region 
III is predominantly composed of the majority spin compo- 
nent. In this regime, one expects the cloud wing-shape to be 
that of a non-interacting Fermi gas, and this provides the ba- 
sis for a reasonable thermometry (JJ. As T is lowered the 



noncondensed pairs in Region II will be converted into the 
condensate, thereby merging Regions I and II. 

Because we have not yet incorporated the qo ^ corre- 
lations of the FFLO state, in Figs, ^andlll the largest of 3 
values of p used is associated with an instability at the very 
edge of the minority cloud. Nevertheless, since essentially 
vanishes there, this is expected to have very little qualitative 
effect on our results. 

The insets in the lower panel show the density profiles 
for the majority (in black) and minority (in red) component. 
Qualitatively similar to what has been observed experimen- 
tally 1 3], a small "kink" in the majority is present at the radius 
at which the condensate ends. The minority component con- 
tains essentially only a condensate central peak with a very 
weak bi-modal structure. The upper panel shows the behavior 
in the normal state, where the condensate A sc — 0. Nev- 
ertheless, it can be seen that an excitation gap A is present 
throughout the cloud. In this way the particle profiles do not 
correspond to those of a non-interacting gas, and the polar- 
ization is rather evenly distributed at all radii in the cloud. It 
may also be seen from these insets that as the system varies 
from above to below T c , the profile of the minority component 
contracts into the center of the trap, as observed @]. 

We present comparable figures for the unitary case in Fig.|2] 
Most of the observations made above for the near-BEC case 
obtain at unitarity as well. Here, one can see from the insets, 
however, that the kink in the majority profile is less apparent. 
Finally, we turn to Fig.|3]where the counterpart plots are pre- 
sented for the BEC regime with 1/kpa = 3.0. From the right 
column, it may be seen that polarization penetrates down to 
the trap center, when the overall polarization p is high. This 
is in agreement with the expectations from the homogeneous 
case 1 11|. 

It should be stressed that our calculations indicate that dra- 
matic changes in the shape of the density profiles do not oc- 
cur until T is substantially lower than T c . This may explain 
why the experimentally observed T c values are less |3| than 
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those we compute. Also important is the fact that the mixed 
normal phase we find here (Region II) is not related to that 
introduced at T = in other theoretical work ifT^ fT9ll . The 
noncondensed pair contribution we consider has no counter- 
part at T = 0. Moreover, there is no evidence from T = 
BdG investigations 1 Lj, 1 1 411 for such a mixed normal phase, 
although whether it is consistent with the FFLO state in a trap 
bears further investigation. 

Essentially all the qualitative observations reported in this 
paper correspond to their counterparts in Ref. 1 3] with one ex- 
ception. In Ref. it is claimed that only regions I and III 
are present in the near-BEC regime, whereas, we find all three 
regions appear, just as in the unitary case. Region II corre- 
sponds quite naturally to the presence of noncondensed pairs, 
expected at finite T in the near-BEC regime. At a more quan- 



titative level, our results at unitarity may change somewhat 
when we include FFLO condensate contributions, and asso- 
ciated higher values of p. While we do not find evidence for 
sharp phase separation as reported in Ref. 1 2], after column in- 
tegration of Figs. 1 1131 a double peaked structure emerges for 
the differenceprofile 5n,2d at low T, as has been claimed ex- 
perimentally jH|2l. Importantly, as stressed in Ref. 01, and, as 
is consistent with our longstanding viewpoint |6], these stud- 
ies show that a significant fraction of the normal region in the 
trap contains strong interactions between the two spin states 
which we associate with the presence of noncondensed pairs 
and related fermionic excitation (pseudo)gap. 
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